Anderson localization or nonlinear waves? A matter of probability 
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In linear disordered systems Anderson localization makes any wave packet stay localized for all 
times. Its fate in nonlinear disordered systems is under intense theoretical debate and experimental 
study. We resolve this dispute showing that at any small but finite nonlinearity (energy) value there 
is a finite probability for Anderson localization to break up and propagating nonlinear waves to take 
over. It increases with nonlinearity (energy) and reaches unity at a certain threshold, determined 
by the initial wave packet size. Moreover, the spreading probability stays finite also in the limit of 
infinite packet size at fixed total energy. These results are generalized to higher dimensions as well. 
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Anderson localization (AL) is a fundamental physical 
wave phenomenon with exponential localization of eigen- 
states in linear wave equations due to disorder (TJ [2] . It 
was originally used to explain metal-insulator transitions 
[3], and has been recently related to experimentally ob- 
served localization of weak intensity light [4J and nonin- 
teracting matter waves [5j on disordered lattice systems. 
The fate of AL in presence of many body interactions 
or corresponding nonlinearity has stayed at frontiers of 
theoretical and experimental physics since 

The pioneering experiments [8j demonstrate both the 
principal possibility to study this effect, but also the se- 
vere length and time scale limitations which restrict a 
quantitative analysis. Theoretical studies become there- 
fore central here. Two faces of Anderson localization are 
usually probed: the absence of the wavepacket spreading 
[TIH2QJ (the focus of this paper) and vanishing conductiv- 
ity lOj. Wave packets, for not too weak nonlinear it ies, 
have been found to spread sub diffusively, disobeying AL 
at least over many time scales fTTHT3] . and remarkable 
progress in studying the regimes and characteristics of 
wave packet spreading has been achieved p!6H2Q] . 

Still, the original questions are open. It remains de- 
bated, whether the observed spreading will continue in- 
finitely or slow down and even stop to restore AL, once 
the wave packet densities become substantially small, as 
conjectured in Ref. [Ml [15] |24] . It is unknown, whether 
there exists a lower bound on the nonlinearity strength, 
beyond which wave packets do obey AL, and therefore 
do not spread at all; or are some time scales diverging 
in this limit? Unavoidable finite size, time, energy, and 
precision limitations will make the most advanced numer- 
ical results not entirely convincing. To achieve a break- 
through, a novel approach is in need. 

The probabilistic description of AL in terms of the 
measure of localized regular trajectories in phase space 
(periodic orbits and tori) appears a promising solution 
for nonlinear disordered systems. Indeed, the mere as- 
sumption that a wave packet is launched on a chaotic 
trajectory leads to the conclusion that chaos remains for- 
ever: Arnol'd conjecture, unproved but widely accepted. 



states the uniqueness of the chaotic region in phase space 
[21]. If initial conditions are chaotic, then the trajectory 
will be unbounded in phase space, corresponding to un- 
limited spreading of a wavepacket. If they belong to a 
localized periodic orbit or torus, then the spreading is 
absent. 

The progress has been quite limited up to now due to 
the difficulty of the task. Persistence of tori has been 
proved for a special class of infinite systems [22j and for 
finite tori dimensionality [23j only. The non-zero proba- 
bility of existence of infinite-dimensional tori in generic 
case (for small enough energies) has been conjectured by 
empirical and numerical arguments in Ref. [24]. Quanti- 
tative estimates for AL probabilities are lacking. 

In this paper we show the probabilistic nature of AL 
in nonlinear disordered systems. At any small but finite 
nonlinearity (energy) value there is a finite probability 
for AL to break up and for nonlinear waves to propa- 
gate. It increases with nonlinearity (energy) and reaches 
unity at a certain threshold, determined by the initial 
wave packet size. Moreover, the spreading probability 
stays finite also in the limit of infinite packet size at fixed 
total energy. This finite value is between zero and one 
for quartic anharmoncity (as realized for two body inter- 
actions and in optical Kerr media), strictly one for the 
lower degree of anharmonicity, and has zero as the lower 
bound for higher degrees of anharmonicity. 

We first address the paradigmatic Frohlich- Spencer- 
Wayne (FSW) type d-dimensional classical lattice [22] . 
for which the linear eigenmodes of a disordered system 
are compact single-site excitations with random frequen- 
cies, and where infinite dimensional Kolmogorov-Arnold- 
Moser (KAM) tori [22] exist. This model can be con- 
sidered as the strong disorder limit of generic classical 
Klein-Gordon (KG) and semi-classical discrete nonlinear 
Schrodinger (DNLS) arrays, where the existence of KAM 
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tori is still debated 



]. Its Hamiltonian reads 



meD{n) 



(1) 

where Xn is the displacement of the n = (ni, . . . , n^^)- 
th particle from its original position, pn its momentum, 
D{n) is the set of its nearest neighbors, and the random 
uncorrelated Sn G [1/2,3/2] are uniformly distributed. 
Unless explicitly specified, a chain d = 1 with quartic 
anharmonicity 7 = 4 is considered. Without the loss of 
generality we assume £0 = 1- 

We consider a wave packet with initial width L and 
determine the probability of L-site localized periodic and 
quasi-periodic solutions of ([T]). We start with periodic 
orbits, which are single-site localized solutions, and derive 
conditions of their destruction (when regular trajectories 
delocalize the corresponding wavepackets are assumed to 
spread). Let us construct an exact time-periodic orbit 
of ([1]) localized at n = 0. We apply perturbation theory 
Xn{t) = XlfcLo ^ri\t) in the small-amplitude limit, taking 
Xn\t) = AnCos{t)do^n as the zero-order approximation 
(note that Sq = 1). In first order we find 



x±i{t) = A±i cost, A±i = 



In higher orders it follows 



X±ri{t) = A±riCOSt, A±ri = 
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For the perturbative solution ([3| to converge the ampli- 
tudes must decay: 
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The probability for this condition to hold at site ±n is 
determined by its random on-site potential e±n and the 
oscillation amplitude of its neighbor: 
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The probability to obtain a localized time-periodic solu- 
tion in the infinite chain V = Yl^=i [P''^^ is bounded 
from above by the probability to have decreasing ampli- 
tudes at least in the first neighbors 
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and from below by the probability to get at least an ex- 
ponentially decaying profile in the infinite system 
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where ShcE <C 1 was used in ^ and where E = Aq/2 is 
the energy of the central site. 

It follows from ^ that it is not possible to construct a 
single-site localized {\An/Ao\ <C 1, Vn 7^ 0) time-periodic 
orbit for ^ > 1/3. For smaller energies, the lower bound 
of spreading probability always remains non-zero and 
scales linearily with the total wavepacket energy back- 
ing the arguments in [24j . 

Next we consider L-site localized solutions to ([T]) 
loosely corresponding to tori. Their existence probabil- 
ity is maximized for sparse packets, when the most ex- 
cited sites are separated by intervals of weakly excited 
ones (due to reduced perturbative corrections, or, equiv- 
alent ly, resonance probabilities, see also [24J). We assume 
such a sparse excitation with L sites with the energy E/L 
per site, separated by at least two non-excited sites, as a 
zero order approximation for the perturbation theory. In 
the first order the problem is reduced to L independent 
single-site problems and, with (|6|, the upper bound for 
the localization probability reads 
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Note, that one can proceed to arbitrary high orders as- 
suming sparser zero-order approximations, ultimately ob- 
taining the lower bound analogue to (|8| by exponentiat- 
ing is a monotonously increasing function of L 
with the limiting localization probability for the infinite- 
size packet 



-6kE 



(9) 



Eqs. (|8| and ^ are central results of this work. Note, 
that a wave packet of size L has 2L phase space vari- 
ables to be defined. Therefore, its relevant phase space 
dimension is 2L. The probabilitiy of obtaining regular 
AL states (or not) is defined through the ratio of the vol- 
ume Vl of all points in this 2I/-dimensional phase space, 
which yield AL localization, to the full available volume 
Vl + ^5 where ^5 is the volume of all points which yield 
spreading: Vl = vi/{vi-\-Vs). We can now conclude, that 
for a wave packet of size L no regular AL states are ex- 
pected if the energy density h = E/L > 1/3. But even 
for < 1 /3 there is always a finite probability to observe 
spreading trajectories. Most interesting is that given a 
fixed total energy E, the probability for AL is approach- 
ing a finite value below one in the limit of infinite packet 
size L ([9|. Therefore even in this limit (of zero energy 
density and infinite packet size) there remains always a 
nonzero probability to spread, i.e. Vs/{vi-\-Vs) 7^ in this 
limit! Remarkably also the derivative dpL/dL ~ is 
vanishing as a power law, and not as an exponential for 
large L. Therefore, at variance to the case of AL, there 
is no new length scale emerging. In particular, already 
the first moment {L) obtained with such a probability 
distribution function diverges. 
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FIG. 1: (Color online) Schematic dependence of the probabil- 
ity Vv for wave packets to stay localized (dark area) together 
with the complementary light area of spreading wave packets 
versus the wave packet volume V (either initial or attained at 
some time t) for three different orders of nonlinearity 7 < 4, 
7 = 4 and 7 > 4. 



We generalize our results to the case of arbitrary de- 
gree of nonlinearity 7 and lattice dimensionality d. We 
consider localized solutions as a sparse wave packet of V 
equally excited sites and characteristic size L and derive 
the upper bound for their existence probability Vy- As 
its volume scales as V ex we straightfowardly arrive 
at the probability to find an AL wave packet 
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We find that for any 7 no regular AL wave packets are 
expected if the energy density h > h^. For smaller h 
there is always a finite probability to launch a spreading 
wave packet. Most interesting is that the fraction of lo- 
calized AL wave packets vi/{vi + Vs) tends to zero in the 
limit F ^ 00 at fixed £^ for 7 < 4, and tends to unity 
for 7 > 4 (note, that this is only the upper bound for the 
AL probability), as shown schematically in Fig [I] 

To probe our analytical results we simulate the evolu- 
tion of single-site initial excitations with energy £^ in ([T]) 
and (11) for (i = 1, 7 = 4. We use a symplectic SABA- 
type scheme p^B) [201 final times t^nd = 10^- To 

characterize the wave packet we calculate its second mo- 
ment m2 = ~ ^^En/E^ with n = ^nEn/E which 

n n 

measures the squared distance between both packet tails. 



and the participation number P 



which tells the 



number of most strongly excited lattice sites. 

We first choose an energy E = 1 which according to 
our above results should yield spreading wave packets 
with probability one. Indeed in Fig|2]we plot the growth 
of 7712 for two disorder realizations, which survives a fit 
m2 cx \^ accord with theoretical preditions p!8l 
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FIG. 2: (Color online) Second moment of a spreading single 
site excitation for two different disorder realizations versus 
time, (i=l,7 = 4, ^ = 1. Dashed and dash-dotted lines 
with the slopes 1/3 and 1/2, respectively, guide an eye. Inset: 
Second moment of a spreading and a nonspreading single site 
excitation for two different disorder realizations versus time, 
(i= 1, 7 = 4, £; = 0.05. 
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FIG. 3: (Color online) Sorted participation numbers Pi for 
Nr — 1000 different disorder realizations 1 < i < Nr at t = 
10^ 5 • 10^ 10^ (bottom-top) for E = 0.08 (bl), E = 0.05 (v), 
E = 0.03 (b), E = 0.02 (g), E = 0.002 (o). Inset: single-site 
localized fractions (Vl^i): numerics by symbols, linear fits 
for FSW and for KG by solid lines. 



19j. For small energy E = 0.05 we expect to observe 
both spreading and nonspreading wave packets. In the 
inset in Figj2]we show two typical cases - in one case the 
second moment does grow, in the second case it does not, 
therefore indicating AL. 

In order to quantify the above results, we perform sim- 
ulations of single site excitations at a given energy E for 
Nr = 1000 different disorder realizations. We measure 
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the participation number P{t) at three large time values 
t = 10^, 5 • 10^, 10^. For each time we sort the participa- 
tion numbers by their value and plot them in Figj3]versus 
the sorting index i. We find that for very small energy 
E = 0.002 all three curves for different times overlap and 
have small values of P, indicating AL at least up to the 
largest time. For larger energies, e.g. £^ = 0.05 we find 
that the curves overlap up to some largest index im{E)^ 
and differ for larger indices. 

To quantify the results we introduce the following nu- 
merical criterion for a single-site (or a single-mode) lo- 
calized trajectory: the one with participation number 
P{t) < 1 + a, choosing a small enough a = 0.2. As 
P ?^ 1 + 2(£;_i + Ei)/Eq for strongly localized tra- 
jectories, ([5| yields a theoretical estimate for the non- 
spreading fraction: V ^ 1 — Q^2/aE. The numerically 
obtained dependence of the non-spreading fraction on the 
energy shows a linear decay and agrees well with this an- 
alytical estimate (Figjsj inset). 

A similar approach can be developed the Klein- Gordon 
(KG) chain described by the Hamiltonian 
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This is a generic model with a finite localization length 
of Anderson modes [H [H HZl [H [20]. We performed a 
similar numerical analysis at = 1 as for the FSW model 
with launching a single eigenmode of the corresponding 
linear eigenvalue problem of u = 0. The parameter W = 
6. The final outcome is plotted in the inset in FigjS) it 
shows again a linear decay of Vl with energy, similar to 
the FSW model. 

In summary, we developed the theory of AL in non- 
linear disordered systems in zero-temperature limit. In 
contrast to linear systems, where AL forces any wave 
packet to stay localized for all times, the fate of AL be- 
comes a matter of probability and nonlinear waves may 
propagate. With increasing the strength of nonlinearity 
(energy) the probability of AL is reduced, and reaches 
zero at a certain level of nonlinearity depending on the 
initial wave packet size. At fixed total energy, the frac- 
tion of AL wave packets does not reach unity in the limit 
of an infinite packet size, but stays finite yet less than one 
for quartic anharmonicity. For weaker degree of nonlin- 
earity the AL fraction in this limit vanishes completely, 
while for stronger nonlinearity degree its upper bound 
tends to one. It follows that the previously conjectured 
slowing down or halt of wave propagation [14 , 15 , 24J are 
not realized at least for the quartic and weaker anhar- 
monicity. The observation of a finite fraction of spreading 



wave packets at small energies questions the correctness 
of the perturbation calculations in [25j, which seem to 
exclude such a possibility. These results generalize to 
higher dimensions as well. An analytic treatment of the 
generic KG-type system yields similar results and will 
be presented in more detail elsewhere. Our results not 
only resolve a fundamental theoretical problem in nonlin- 
ear wave physics but also have multifaceted experimental 
implications, as for design, as for interpretation. 

MI acknowledges financial support of the Dynasty 
Foundation. 
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